QCD at non-zero temperature and density from the Lattice 
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The study of systems as diverse as the cores of neutron stars and heavy-ion collision experiments requires the 
understanding of the phase structure of QCD at non-zero temperature, T, and chemical potential, fi q . We review 
some of the difficulties of performing lattice simulations of QCD with /j, q 7^ 0, and outline the re- weighting method 
used to overcome this problem. This method is used to determine the critical endpoint of QCD in the (/j, q ,T) 
plane. We study the pressure and quark number susceptibility at small fi q . 



1. Introduction 

The QCD phase diagram is a rich source of in- 
teresting phenomena which is key to understand- 
ing a wide variety of physical systems, e.g. the 
early universe, the cores of neutron stars and 
heavy- ion collisions experiments. Figure ^ illus- 
trates the QCD phase diagram in the chemical 
potential (fj, q ) temperature (T) plane. It depicts a 
very rich structure including familiar phases such 
as nuclear matter, together with very unfamiliar 
ones such as the colour superconducting phase. 
See e.g. £Q for further information. 

Lattice simulations of QCD with non-zero 
chemical potential have recently provided impor- 
tant qualitative and quantitative information on 
the QCD phase diagram in the (T,/j, q ) plane 
This work will concentrate on studies at finite 
T and small \x q which are directly applicable to 
heavy-ion collision experiments. These are shown 
as crosses in fig Q] labelled by the name of the rel- 
evant experiment. 

We begin this paper by outlining the difficulties 
of simulating QCD with non-zero /i q , and then de- 
fine the reweighting method which is the basis for 
much of this work. The method is then used to 
determine the crossover transition near /i = be- 
tween the confining and deconfining phase (which 
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Figure 1. The QCD phase diagram in the 
temperature-chemical potential plane. First or- 
der phase transitions are shown as solid lines, and 
second order transitions by filled circles. Crosses 
depict heavy-ion collision experiments. The re- 
gion of interest in this work is the crossover tran- 
sition at T r w 170 MeV. 



is depicted as a filled circle in fig^ncar the RHIC 
point). We also use this method to determine 
the pressure and quark number susceptibility for 
QCD. Full details of this work appear in [3] & 0]. 
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2. The Sign Problem 

In lattice simulations of QCD at \i q = 0, the 
quark matrix, A4(fi q = 0) = Ip+m, has the prop- 
erty 



M(p q = 0)= 75^0^ = 0)75 



(1) 



which means that the Boltzmann weight is real, 
and can be interpreted as a probability measure. 
However, when a non-zero chemical potential is 
introduced, the quark matrix becomes (schemat- 
ically) M(fi q ) — ]p+ m + fi q jo- K therefore no 
longer has the hermiticity property, and the cor- 
responding Boltzmann weight is complex (the so- 
called Sign Problem), removing the probabilistic 
interpretation. This means that the traditional 
Monte Carlo approach of estimating the path in- 
tegral is a non-starter. 

A number of ways have been promoted to over- 
come this problem. These fall into two methods. 
The first is to study models of QCD which, due to 
their simpler mathematical formulation, do not 
suffer from the sign problem for fi q ^ 0. Such 
models include the NJL model and 2-colour QCD. 
(See [H] for more details.) 

The second approach effectively studies QCD 
at n q = + e by using either a reweighting tech- 
nique 6 7 8 3 9 or an analytical continuation ap- 
proach |10lllj . The latter approach simulates 
with an imaginary fj, q and analytically continues 
back to the real fi q — axis. The reweighting tech- 
nique calculates the expectation value of some op- 
erator, (O), by factoring off the \i q ^ part of 
the Boltzmann weight and appending it to the ob- 
servable O. Using this approach ensembles gen- 
erated at jjbq = can be used. Expressing this 
idea mathematically we have, 
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In these equations, S g is the gauge action, U de- 
notes the gauge degrees of freedom, (3 = 6/g 2 , 
where g is the gauge coupling, and we have in- 
troduced the dimensionless lattice chemical po- 
tential, /i = fi q a. By setting /io to zero, eq(|2J) 
enables the calculation of quantities at /i ^ us- 
ing ensembles generated at /iq = 0. 

The results presented in this work are derived 
from the reweighting approach as well as a Taylor 
expansion of {O) about [i q = 0. 

3. Lattice details 

Our simulations were performed on a N^ x 
Nt — 16 3 x 4 lattice with two (degenerate) quark 
flavours (see |3] @] for full details). We used two 
different quark masses, ma = 0.1 and 0.2 corre- 
sponding to pseudoscalar-vector meson mass ra- 
tios of Mps/My ~ 0.70 and 0.85 respectively 
|12j . A Symanzik improved gauge action together 
with the p4 improved quark action ^3] was used 
to reduce lattice discretisation artefacts. 

A range of (3 values was simulated, each corre- 
sponding to a different temperature via the usual 
relationship, T = l/(a(f3)N t ), where Nt is the 
number of lattice sites in the temporal direction. 
0(1, 000, 000) trajectories were generated in this 
project using the APEmille computers in Biele- 
feld and Swansea. 

4. Reweighting preliminaries 

As an illustration of the reweighting method we 
show the susceptibility of the chiral condensate, 
Xz^ as a function of (3 in fig [21 for a quark mass 
of ma = 0.1. Note that to produce this plot, a 
reweighting has been performed both in f3 and fi 
(i.e. the simulations were only performed a rela- 
tively coarse values of (3 compared to those shown 
in fig©. 

As can be seen, there is a clear peak in the 
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Fi gure 2. The chiral susceptibility, x^"^ 
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susceptibility, which occurs at the crossover tran- 
sition. This means that the expectation value 
of tjnj) fluctuates strongly at this (3 value corre- 
sponding to a temperature T = l/(a(f3)N t ). By 
mapping out the position of this transition as a 
function of f3 and n, the phase diagram, figOl can 
be obtained 3 . 

Note that figs 01 and were obtained by Tay- 
lor expanding the susceptibilities to 0([i 2 ). The 
phase diagram, fig [3] is therefore correct to this 
order in fi 

5. Pressure 

The pressure is related to the grand partition 
function, Z(T, V,fi q ) by 



P_ 
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VT 3 



InZ. 
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We can now Taylor expand both sides of this 
equation in the dimensionless ratio (Aq/T. Since 
we can only utilize Monte Carlo methods at /i q = 
0, we perform this expansion about this point ob- 



Figure 3. The phase diagram obtained from our 
reweighting technique. The errors shown are sta- 
tistical. The diamond is the endpoint of the first 
order phase transition obtained by Fodor and 
Katz [7]. 



taining 
(JLi 



) = - 




p_ 




J - T4 




J 1 ! 





1 [i\ 9 2 (p/T 4 ) 1 $ 9 4 (p/T 4 ) 



2! T 2 d(fi q /T) 2 

oo 



4! T 4 d(fi q /TY 



n=2,4, 



(4) 



where all derivatives are taken at fi q — 0. Only 
even powers occur in eq(@J because odd deriva- 
tives of the free energy with respect to fi q vanish 
U3- Using eq©, the derivatives in eql@} can be 
expressed as expectation values (over a fi q = 
ensemble) . 

In figs 0] and [5] the first two coefficients of p, 
C2,4, are plotted as a function of temperature. 
For comparison, the Stefan-Boltzmann (SB) val- 
ues are plotted. These are valid in the limit of 
infinite temperature, where the system becomes 
asymptotically free. The continuum SB expres- 
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sion for massless quarks is 
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so the coefficients 02,4 can be simply read-off. 
Also shown in figs 0] and are SB values cor- 
responding to a free field lattice calculation at 
Nt — 4 (i.e. corresponding to the temporal ex- 
tent of our lattice, see sec. |HJ). The details of this 
calculation can be found in 4 . As can be seen 
from figs 0| and (Sj both C2 and C4 vary sharply in 
the critical region, T sa Tq. 
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Figure 4. The quadratic coefficient of the Taylor 
expansion of the pressure, C2, as a function of 
temperature. (See eq@.) 



Using the coefficients €2,4, we plot the pressure 
difference, A(p/T 4 ) as a function of temperature 
for various n q /T in fig El Again we observe an 
abrupt change in the pressure as the critical re- 
gion is crossed. It is interesting to note that the 
pressure at finite density is very similar to its 
fi q — value. For instance, at the RHIC point 
(fiq/T m 0.1) the pressure is only 1% more than 
its zero density value. 
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Figure 5. The quartic coefficient of the Taylor 
expansion of the pressure, C4, as a function of 
temperature. (See eq(@J.) 



6. Susceptibility 

We turn our discussion to the quark number 
density, n q , and its susceptibility, \qi defined as, 



Xq 

J>2 



d 2 ( P /T 4 ) 



The ingredients required to calculate the suscep- 
tibility have already been determined in the Tay- 
lor expansion of the pressure (see eqs(Pl & EJ^- 
Figure plots Xq/T 2 as a function of tempera- 
ture for various fi q /T values. As can be seen, a 
peak in \ q develops at T w To which grows with 
fi q . This is a signal of a critical endpoint in the 
(p q ,T) plane where the fluctuations in the num- 
ber operator ipjoif) grow. Also note that the peak 
occurs at T = T Q , implying that the chiral sym- 
metry restoration transition observed in fig[21and 
the deconfinement transition of fig [7] coincide P^. 
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Figure 6. The dimensionless pressure difference 
A(p/T ) as a function of temperature for various 
Hq/T values. 



7. Conclusions 

We have presented some results of studies of 
QCD at non-zero temperature and small chemical 
potential j-'il 1 j . In particular we have summarised 
the difficulties in using conventional Monte Carlo 
techniques to perform lattice simulations of QCD 
at non-zero [i. We outlined the reweighting tech- 
nique which has recently been very successful in 
obtaining numerical results at fj, ^ 0, despite this 
difficulty |7I3| . This method was used to deter- 
mine the position of the crossover transition be- 
tween the confined and de-confined phases. Fi- 
nally, results for the pressure and susceptibility 
were produced, based on a Taylor expansion in 
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